Components of multifractality in the Central England Temperature anomaly 
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We study the multifractal nature of the Central England Temperature (CET) anomaly, a time series that spans 
more than 200 years. The series is analyzed as a complete data set and considering a sliding window of 11 years. 
In both cases, we quantify the broadness of the multifractal spectrum as well as its components defined by the 
deviations from the Gaussian distribution and the influence of the dependence between measurements. The 
results show that the chief contribution to the multifractal structure comes from the dynamical dependencies, 
mainly the weak ones, followed by a residual contribution of the deviations from Gaussianity. However, 
using the sliding window, we verify that the spikes in the non-Gaussian contribution occur at very close dates 
associated with climate changes determined in previous works by component analysis methods. Moreover, the 
strong non-Gaussian contribution found in the multifractal measures from the 1960s onwards is in agreement 
with global results very recently proposed in the literature. 
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I. INTRODUCTION 

After Mandelbrot setting forth the concept of fractal, 
the interest of the scientific community in systems ex- 
hibiting self-similarityi, i.e., the property by which a 
system is (approximately) equal to a part of itself has 
soared. Since then, fractals have been applied to such di- 
verse fields as physiology and economics with the inten- 
tion of analyzing the existence of scale-invariant behav- 
ior. This property can be characterized by a sole fractal 
(or Haussdorf ) dimension or a spectrum of locally depen- 
dent exponents, the so-called multifractals spectrum^—. 
Considering a given observable O, the scale invariance 
(self-similarity) we have mentioned can be mathemati- 
cally written as, 



f({\O z }) = \ a ^ f({0}). 



(1) 



Time series are generally self-affine, i.e, they have differ- 
ent scaling properties in in x and y directions. When, 
and only when, the exponent a(z) is constant for all z, 
the system has single scale invariance behavior. 

In the field of meteorology and geophysics, fractality 
and scale invariance are popular concepts as well. How- 
ever, the large majority of the surveys over self-similarity 
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of climate time series have been only devoted to the ap- 
praisal of memory effects^— introduced in the form of 
linear correlations or correlations that can be written in 
a linear way. The existence of memory is quantified by 
the power spectrum exponent /3, which is related to just 
one specific value a (the so-called Hurst exponent) of the 
whole multifractal spectrum. This single a is related to 
the power spectrum exponent (3 via the Wiener-Khinchin 
relation and to the fractal dimension, D, by D = 2 — H . 
When a(z) ^ constant, it is possible to extend these 
studies beyond standard memory effects in order to ex- 
plore and perceive a broader range of statistical and dy- 
namical properties. 

In this manuscript, we shed light on the multi-self- 
similar nature of a paradigmatic meteorological time se- 
ries: the anomaly of the Central England Temperature 
(CET), which is the largest instrumental daily data set 
available on the subject. Our work aims at providing re- 
liable quantitative answers to the following questions: i) 
To what extent does the CET anomaly series evince real 
multifractality? ii) What is the weight of each contribut- 
ing mechanism to the measured self-invariant signature? 
Hi) What is the evolvement of those contributions and 
what are the possible pivotal changes in the nature of 
the system? 

The present manuscript is organized as follows: in Sec. 
77 we provide a thorough description of the data analyzed 
and the methods we employed; in Sec. Ill we present our 
results concerning the proposed questions, and in Sec. 
IV we discuss the conclusions we have reached from our 
enquiry. 
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II. DATA AND METHODS 

A. The Central England Temperature time series 

The data of the Central England Temperature (CET) 
is the longest instrumental time series, starting in 1772. 
Owing to its large span, it can be used to perform com- 
parisons between periods in which greenhouse gases emis- 
sions were lower and contemporary times. The CET cor- 
responds to a daily average of measured temperatures at 
different climate stations close to the triangle defined by 
the English towns of Bristol, Lancashire and London (see 
Fig. 
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FIG. 1. Excerpt of Great Britain's map with climate stations 
locations according to the legend shown in the figure. © 
Crown 2009, the Met Office (Reproduced with permission.) 

The first compilation of a monthly series was presented 
by Manle y 10 ' 11 covering the years from 1659 to 1973. Af- 
terwards, these data were updated to 1991 when Parker 
et at (1992) calculated the daily series. Paying heed 
to urban warming, the data have been adjusted by 0.1- 
0.3 Celsius degrees since 1974. From January 1878 un- 
til November 2004, the data was recorded by the Met 
Office using the stations of Rothamsted, Peershore Col- 
lege and Ringway and from the latter date onwards the 
station of Stonyhurst replaced Ringway. Concomitantly, 
revised urban warming and bias adjustments have been 
applied to the Stonyhurst data. In our work, we specif- 
ically set our sights on analyzing the difference between 
the daily temperature and the average temperature cycle 
temperature obtained by determining the average over 
all filed years of the temperature for each day making 
up a year in the Julian calendar. The data used in this 



work have been obtained from KNMI (The Royal Nether- 
lands Meteorological Institute) Climate Explorer website 
http : //climexp . knmi . nl/ where besides the daily tem- 
perature, the temperature anomalies and the seasonal cy- 
cles are also provided. The annual average of these time 
series is presented in Fig. [5J 
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FIG. 2. (Color online) Temperature anomaly (annual aver- 
age) for the Central England Temperature (CET) time series. 

We proceed twofold; First, we analyze the overall self- 
similar characteristics of the CET anomaly and after- 
wards we pore over the non-stationarities of the time se- 
ries. Specifically, we scan the daily time series, from 1772 
to 2007, with a 11-year sliding window, in which we ap- 
ply the same methods that we describe hereinafter. Each 
window starts on the 1st January of year i and ends on 
the 31st December of the year i + 10, resulting in window 
of 11 x 365 points. Hence, there will be N— 10 (226) win- 
dows whose centers are located at the year i + 5, where 
N is the number of years in the series. We have cho- 
sen this sliding window length for the following reasons: 
a) It has approximately 4fc data points, which is enough 
for an accurate estimation of multifractal spectra; b) It 
produces a sufficient number of windows for an appropri- 
ate evaluation of the time evolution of scaling properties 
of the system, since using of larger windows reduce the 
number of them; c) It does not affect our spectra and 
thus our results, since no spurious scale with frequency 
of 1/llyears was found. In other words, using over- 
lapped instead of simply juxtaposed windows, any sig- 
nificant change that might miss in one windows will be 
detectable in the next one; d) The length of 11 years is 
large enough to cope with El Nino (La Nina) phenomena. 

B. Methods 

1. Multifractal detrended fluctuation analysis - MF-DFA 

Detrended fluctuation analysis (DFA) has been 
pointed out as a tool to monitor climate^ and such 
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methodology has been employed in the study of the time 
evolution of the Hurst exponent of temperature time se- 
ries, by using of moving windows 15 . Briefly, it corre- 
sponds to a method which permits quantifying the way 
the variance of the cumulative signal scales with the 
size of the time interval considered for aggregation, af- 
ter a detrending procedure that seeks to remove non- 
stationarities of the signal. This method equally weighs 
large and small deviations and is thus unable to fur- 
nish any information on the distinct influence of both 
cases. Notwithstanding, this method can be easily mod- 
ified to determine unequal influence of large and small 
variations or, in other words, to appraise multi-scaling 
properties: the multifractal detrended fluctuation analy- 
sis, MF-DFA. 

The MF-DFA 16 is one of the most applied methods 
to determine the multifractal properties of time series 
in several field o 17 ' 18 . We have chosen to apply the MF- 
DFA rather than the Wavelet Transform Modulus Max- 
ima (WTMM) 19 taking into account a recent compara- 
tive study in which it has been shown that, in the ma- 
jority of the situations, the MF-DFA presents reliable 
results 20 . Specifically, the MF-DFA does not introduce 
as specious multifractality as the WTMM. 

Considering the time series {x (t)} (x represents the 
temperature fluctuation in our case) made up of N ele- 
ments, the MF-DFA is composed of the following steps: 

• Determine the profile Yi (t) that corresponds to the 
deviation of the signal elements from the mean 

t 

Y t (t) = (*) ~ (*>] . (l<t<N), (2) 

1=1 

where (...) represents the average over elements 
and i = 1, N. 

• Divide the new profile Yi (t) into N s = int (— ) non- 
overlapping intervals of equal size s; 

• Compute the local trend by a least-squares adjust- 
ment method, and thereupon the variance, 

F 2 ("> s) = \Y j {Y[{ V -l)s + l]- y v , (3) 
i=i 

for each segment v = 1, . . . , N s , where y v {%) repre- 
sents a p*' l -ordcr fitting polynomial in the segment 
v. Since the number of points in most time series is 
not a multiple of the number of segments in a given 
scale s, the same procedure may be performed on 
the series with the order of its elements reversed. In 
this case, the two resulting matrices must be con- 
catenated and the fluctuation function described in 
the next item must be divided by two. The order of 
the polynomial is relevant for the results one might 
obtain. For the CET anomaly time series we used 
polynomials of 2 nd -order. 



Determine the average F z (s) over all segments to 
obtain the fluctuation function of order z, 



l/z 



and 



^(s)^cxp j^L^ln [p 2 (,/,«)] j 



V z#0 , (4) 



0. (5) 



Equation (jl]) is basically the only equation different 
to the DFA method. On the one hand, negative 
values of z decrease the influence of large values 
of F 2 , whereas it augments the influence of small 
values. On the other hand, positive values of z have 
exactly the opposite effect. It is from this outcome 
that we are able to check the different scales present 
in the system. 

• Assess the scaling behavior of F z (s) considering a 
log — log scale representation of F z (s) versus s for 
each value of z. If the series {x(t)} shows multi- 
scaling features then, 



(6) 



When the exponent h(z) is negative, as in the case 
of anti-correlated signals, or close to zero, the double 
summation 



(7) 



1=1 



must be used in Eq. ([3]) in lieu of Eq. ([2]). This ad- 
ditional step is necessary because the methodology as it 
was presented is inaccurate for negative generalized Hurst 
exponents. The generalized fluctuation function (Eq. [6]) 
becomes then 



Fz (s) 



(8) 



where h(z) — h(z) + 1. 

Small fluctuations are generally characterized by large 
scale values of the exponent h (z) (and z < 0), whereas 
large fluctuations are typified by small values of h (z) 
(and z > 0). 

Bridging this procedure and the standard formalism, 
we verify that [F z (s)] z can be interpreted as the partition 
function, Z z (s)pl which is known to scale with the size 
of the interval as, 



Z z (s) 



r(z) 



(9) 



Hence, according to Eq. and Eq. (O we have, 

r(z) = zh{z)-l. (10) 
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Legendre transforming, 

f(a) = za-T(z), (11) 

we can relate the exponent r (z) with the exponent, a, 

, , . dh(z) , „. 

a = h(z) + z—±-t, (12) 

and 

f(a)=z[a-h(z)] + l. (13) 

For z — 2, h (2) = iJ corresponds to the Hurst expo- 
nent customarily determined by methods like the original 
R/S ratio or the DFA 2 ^ w ith H = (/3 + l)/2 and /3 is 
the power spectrum exponent according to the Wiener- 
Khinchin relation. On the other hand, z = give us 
the support dimension of the self-afHne structure. In the 
case of a monofractal, h (z) is independent of z implying 
homogeneity in the scaling behavior and thus Eq. (TTU)) 
reduces tor(z) = zH — 1. Explicitly, there exist only dif- 
ferent values of h (z) if large and small fluctuations scale 
in differently. Lastly, we must stress out we are using the 
term fractal in not totally formal way. In fact, despite 
the fact that in some cases h{z) to the box-counting di- 
mension, it generally does not verify all of the properties 
to be classified as the Hausdorff-Bcsicovitch dimension of 
the signal. 

2. Components of multifractality 

There are two ingredients that induce multifractality 
in a time series: the deviations from a scale dependent 
distribution and the dependence between its elements^. 
These factors are generally assumed as independent and 
thence the total multifractality usually stems from the 
superposition of both. This means we can accredit a cer- 
tain value /ing( z ) to the non-Gaussianity and a value 
h(z) — /ing(z) to the dependence. Q However, it should 
be noted that this is not necessarily true. For instance, 
in the case of heteroscedastic models such as ARCH-likc 
proposals the two elements of multifractality are strongly 
dependent^. This assertion is based on the fact that the 
non-Gaussianity of the time series comes about from the 
existence of dependencies between elements that are lin- 
ear correlations in the (square) volatility. When these 
dependencies are set to zero, the outcome is a Gaussian 
series. Consistently, taking into consideration that the 
variables are uncorrelated, we verify that the two contri- 
butions towards multifractality are one and the same. 
From Eqs. (TIl3|) - (fl"3)) . we can understand that a plain 



Throughout the text we employ the term non-Gaussianity some- 
what inaccurately to refer to a distribution which does not have 
a defined scale as the Gaussian presents. Nevertheless, it is a 
more or less established terminology for the absence of scale. 



way of checking the existence of multifractal features is 
achieved by computing the difference, Aa, between the 
minimum and the maximum values of the spectrum f(a). 
If there is a mono-fractal nature we shall have Aa = 0, 
otherwise we should be able to separate out the contri- 
butions of the non-Gaussian deviation, A«ng, the non- 
linear dependencies, Achnl, and linear features, Aculd, 
so that 

Aa = Aong + Aculd + Ao;nl- (14) 

It is generally assumed that linearities do not con- 
tribute to the multifractal broadness because they are 
traditionally related to a single scale of dependence. 
However, previous results on financial dat a 18 ' 25 have 
showed that this linear/weak dependencies do modify the 
MF-DFA spectrum and thus we shall a priori not scrap 
its contribution. Consider that our time series {x (t)} 
is associated with a non-Gaussian distribution and that 
its elements exhibit dependencies at a non-null level. If 
we perform an appropriate shuffling process, we define 
a new series related to original one by having the same 
probability density function, but for which there is no 
dependency between the elements, because the memory 
was totally destroyed. As proven in ref. 1 161 when we 
analyze the multifractal spectrum of the shuffled surro- 
gate, the curve will depart from the mono-fractal behav- 
ior, Aa = Aa( s huf) 7^ 0, due to the non-Gaussianity of 
the probability density function and thus, 

Aa NG = Aaff huf) , (15) 

where Aa cff is an effective Aa that will be defined 
shortly. Alternatively, we can apply a procedure of phase 
randomization in which we Fourier transform the series in 
the uj domain [— 7r, tt) and replace the phases by random 
values (uniformly distributed) up to the first half of the 
transformed series (excluding the phases of to = — tt, 0). 
The conjugates of the first half are then used for the re- 
maining terms. Afterwards, we apply the Inverse Fourier 
transform. The final result is a time series in which the 
previously existing deviations from the Gaussian distri- 
bution are destroyed (see an example in Fig. [3J. Nonethe- 
less, we can verify that the power spectrum of the sur- 
rogate is the same as the original series. Accordingly, 
since the power spectrum describes the existing linear- 
correlations, when we evaluate the multifractality of the 
new series, Aa( ran( j), the only contribution originates 
from the linearities of the system, 

Aa LD =Aa£ nd) . (16) 

When we analyze the multifractal nature of a time se- 
ries which has been shuffled and phase randomized (or 
vice-versa) we are supposed to get Aa( s huf&rand) = 0. 
Nevertheless, the values different from zero (consider- 
ing the error bars) that are found when we compute 
Aa( s huf&rand) are still important as they can be used to 
check the effect of the error introduced by the finiteness 
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FIG. 3. (Color online) The (□) symbols were obtained from 
a series of 10 4 elements following a (t = 3)-Student distri- 
bution (or a (q = 3/2)-Gaussian with a 2 = 3) which was 
shuffled in order to remove any possible dependence between 
the variables due to the pseudo-random generator and the 
(•) symbols were obtained from a surrogate time series orig- 
inate sprung from the phase randomization procedure which 
matches a Gaussian distribution with the same standard de- 
viation. 



of the time series and the method, 
define the effective Aa as: 



This allows us to 



Aaf, 



Aa 



(...) 



Act 



(shuf&rand) i 



(17) 



for the original, shuffled and random phase time series. 

We can make an appraisal of the joint contributions 
of the non-Gaussianity and the non-linearities as it be- 
tides ARCH-like processes by proceeding the following 
way. Instead of randomizing the phases of the Fourier 
Transform x(lj), we preserve them and assign a con- 
stant value to the absolute value of x(u>) heeding that 
\x(u>)\ = \x(—w)\. With this, we are setting a constant 
power spectrum typical of a white noise while preserv- 
ing every other feature of the dynamics (by preserving 
the phase). After Inverse Fourier transforming, we can 
compute the multifractal broadness Aapp, the effective 
value of which is to be compared with the sum of Acung 
and Aonl- 

In spite of the fact that the effect of time dependencies 
in the multifractality is already take into consideration 
in the MF-DFA procedure (The D stands for detrended), 
it is worth mentioning that in order to apply the Wiener- 
Khinchin relation the series must be stationary or at least 
close to it. With the purpose to overcome the problem 
of the non-stationarity of the time series for a proper 
evaluation of S(u), which is only necessary to generate 
the surrogate series, we have applied a high-pass filter — 
which removes the non-stationarities. We have also used 
the Burg algorithm 28 i 29 for estimating S(ui) since it is 
very difficult to distinguish noise from information from 
the standard FFT spectrum. 



III. RESULTS 
A. Overall results 

In this sub-section we first analyze the time series as 
a whole set. The results of the multifractal analysis of 
the CET anomaly time series and its various surrogates 
described in the previous section are depicted in Fig. |U 
For all the curves we verify the fat-fractal nature of the 
series as the maximum value of /(a) is equal to 1. By ex- 
trapolating the values a m ; n and a max at which the /(a) 
curve intersects the a axis, we determined the values of 
Aa = 0.37 for the original times series and for the shuf- 
fled plus randomized surrogate the multifractal broad- 
ness Aa( s h u f& r and) = 0-16 yielding an effective broadness 
Aa off = 0.21. Concerning the remaining surrogates we 
have obtained effective values Aa^ uf j = AaNG = 0.02, 
^ a (rand) = ^ a LD = 0.13 and AaNL = 0.06. Examin- 



ing the multifractal results of the phase-preserved sur- 
rogate we obtained Aa( w m te ) = 0.05. Taking into ac- 
count the error (±0.01 in our analyses) we verify that 
Aapp ~ Aq;nl + Aoing- 
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FIG. 4. (Color online) Multifractal spectrum f(a) vs a for 
the original CET anomaly and for the surrogates generated 
by shuffling, phase randomizing, shuffling plus phase random- 
izing and phase preservation described in Sec. Ill B 21 . 

The slender contribution of the Non-Gaussianity is eas- 
ily understandable by the well-known nearness to the 
Gaussian distribution of the CET anomaly. 

These results, particularly the influence of the linear 
correlations can be double-checked the other way round. 
Specifically, we can create a new surrogate which retains 
the power spectrum S (to) (shown in Fig. [3]) of the CET 
anomaly but whose elements are associated with a Nor- 
mal distribution and which do not present any non-linear 
dependencies. This can be made by generating a Normal 
distributed and independent time series that is ulteri- 
orly Fourier transformed. The amplitude of each ele- 
ment of the set {x (w)} is changed by multiplying it by 
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FIG. 5. (Color online) Power spectrum S(ui) vs uj for the CET 
anomaly obtained using the Burg algorithmi-^i-^ The dashed 
lines represent the error margins ( 5 % of the power spectrum 
value) . 



W S (lu). Thereafter, we apply the inverse Fourier trans- 
form on this series. The results of this procedure are 
depicted in Fig. which shows an effective broadness 
^ a (fouricr) = 0-15, completely compatible with the result 
Aculd hereinabove mentioned. 
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FIG. 6. (Color online) Multifractal spectrum f(a) vs a for 
the original CET anomaly and for the surrogates in which 
the power spectrum is kept. It is visible that the broadness 
of both spectra are not far away. Moreover the broadness 
Aa (Fouricr ) falls within Aa (rand) . 



features in the series. The first of all indicates the non- 
stationary nature of the multifractal spectrum in the time 
series as depicted in Fig. [71 This multifractal evolvement 
can be tweaked in order to obtain the effective multi- 
fractal broadness Aa eS . The results of this tuning are 
striking. In average we verified that the effective multi- 
fractality represents just about 50% of the broadness ob- 
tained in direct measurements of the CET anomaly time 
series in Fig. 7. Concerning the truthful elements of 
multifractality we can use Eq. (fT5"]) and Eq. (|16D . For the 
contribution of the non-Gaussianity, Aa^G, we under- 
stood that only 33 out of the 226 Aa^^ are actually 
different to zero and consequently to a single-structure 
nature. This corresponds to an average contribution of 
merely 3% which once again is in accordance with the al- 
most Gaussian nature of the time series. Regarding the 
contribution Aa^u, it is obtained when we convert our 
series into a Gaussian series by means of the phase ran- 
domization procedure as we depict in Fig. [3] The aver- 
age contribution towards the multifractal nature is 59%. 
Summing the two contributions we grasp that there are 
missing elements. This implies that the multifractality 
introduced by the non-linearities corresponds to an aver- 
age of 38%, as shown in Fig. [S] 
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FIG. 7. (Color online) Multifractal broadness Aq of the Cen- 
tral England temperature anomaly between 1777 and 2002, 
coarse and effective. The Aa e ^ is obtained after removing 
finite size effects and systematic algorithmic error, by sub- 
tracting from the original Aq, the value of Aa( shll ^&, ran[J ) as 
described in Eq. (|17p . Each curve follows the legend in the 
figure. 



B. Time dependence of the components of multifractality 

Applying the methods described in the previous sec- 
tion on the CET anomaly time series we identified several 



C. Detecting changes in climate 

The results we have presented in the previous sub- 
section can be interestingly compared with the results 
by Berkes et al=. Based on Functional Data Analy- 
sis (FDA), Berkes et al. pointed out a set of climate 
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FIG. 8. (Color online) Non-Gaussian and linear depen- 
dence (correlations) contributions to the effective multifractal 
broadness Aa eff of the Central England temperature anomaly 
between 1777 and 2002. Each curve follows the legend in the 
figure. The Non-Gaussian character is related to the broad- 
ness of the multifractal spectrum of the shuffled time series, 
while the linear dependence is related to the broadness the of 
the multifractal spectrum of the randomized time series. All 
the calculations have been performed using the effective Aa 
given by Eq. l[T7|. 



changes whose major statistical significance occurred in 
1780, 1815, 1926 and 2007 for a previous method called 
MDA and 1780, 1808, 1850, 1926, 1992 and 2007 for their 
new methodology. Looking to our results we can ver- 
ify changes in the non-Gaussianity around 1800 followed 
by flurried periods between 1820-1829, 1834-1845, 1891- 
1897 and a last and standing ruffled period from 1966 
on. Despite each study surveys different quantities, it is 
well-established that changes in a certain dynamics can 
be identified from variations in different observables that 
are associated with that dynamics. It is thus reliable 
that this fair concurrence provides us with an indication 
about the robustness of past changes in the leading dy- 
namical mechanism of the temperature. Although we 
could not observe any change in 1926 or thereabouts in 
Fig. [5] we noticed that there is a sharp peak close to 
this date in Fig. [9l It is worthy to be emphasized that 
the authors of the FDA work refer to the method as a 
"mere modeling assumption that is useful in identifying 
patterns of change in mean temperature curves" . The 
broadest and highest contributions of non-linear depen- 
dence and non-Gaussianity start from the mid 1960s to 
the late 1990s. A particular pattern composed of suc- 
cessive positive values in the annual average of the tem- 
perature anomaly is reported by Parker et al. (1992) 
in this period, which encloses the change of phase that 
is well-documented in the climatological temperature 31 . 
Moreover, we have observed that the spikes verified in 
the A«ng concur with the emergence of kurtosis excess 
when statistical moments are surveyed 3 ^. Recent results 



FIG. 9. (Color online) Non-linear dependence contribution 
to the effective multifractal broadness Aa cff of the Central 
England temperature anomaly between 1777 and 2002. The 
weight of the non-linear dependence contribution are measure 
from the broadness of the /(a) curve of the surrogate time 
series with the same Fourier spectrum of the original one, but 
without linear correlations. All the calculations have been 
performed using the effective Aa given by Eq. {TTJl- 



published by Hansen et al££ using a standard statistical 
approach on a global scale study has verified that the 
temperature anomaly in the last decades are compatible 
with non-Gaussian distributions. Motivated out by the 
time dependence of the scaling properties of the CET 
anomaly we have made an analysis of the a m i n and a max 
(effective values) that we present in Fig. [10] We took note 
of the maxima in the amplitude of both spectra that con- 
form with decadal oscillations. Namely, we have found 
for a m - ln maxima at 16 ~ 1 and 32 " x year" 1 and for a max 
maxima at 17 , 34 year" 1 with the former having an 
extra peak at 55 _1 year -1 . These decadal oscillations are 
not present when the broadness Aa is analyzed. One rea- 
son for this stems from the fact that both a max and a m i n 
show very close peaks in the spectrum and thus when the 
difference between them is considered those peaks are not 
perceived as we have two co-evolving quantities. 



IV. FINAL REMARKS AND OUTLOOK 

In this manuscript wc studied the multifractal features 
of the CET anomaly. This was carried out by consider- 
ing the whole time series and by assuming sliding blocks 
of 11 years. Our results showed that the main compo- 
nent of the multifractal structure of the time series is due 
to the linear dependencies in the dynamics, followed by 
non-linear dependencies and with a residual contribution 
of the deviations to the Gaussianity. Nonetheless, when 
we analyzed the evolvemcnt of multifractality, we veri- 
fied that until 1950 the non-Gaussian contribution comes 
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FIG. 10. (Color online) Amplitude of Fourier transform 
(a(w)) of the maximum and the minimum values of a(t), 
against frequency. Here t represents the center of the slide 
window. The circles identify the local maxima described in 
the text. 



about in periods close to dates at which occurred climate 
changes according to the results obtained by Functional 
Data Analysis^. From the 1960s onwards the contribu- 
tion of the non-Gaussianity becomes significant most of 
the time, a result that matches a recent finding of tem- 
perature anomalies significantly beyond the 3 a criterion 
that characterize the validity (or not) of the Gaussian 
distribution^. 

Finally, although the scope of the present work con- 
cerns the quantitative description of the temperature 
anomaly multifractality and its ingredients as well as its 
connection with climatic changes, our analysis can have 
direct implications in temperature anomaly modeling. 
Our results point out that the temperature variability 
modeling must take into account the multifractal nature 
of the temperature time series, and the contribution of 
each ingredient for the multifractality, as well. A tenta- 
tive dynamical scenario is that wherein we resort to the 
statistical mechanics relation between temperature and 
standard deviation, i.e., the temperature is proportional 
to the standard deviation, and consider a model inspired 
in cascade models of the latter quan tity for turbulent 
fluids such as those introduced in refs. I34l - l36l . Instead of 
considering a cascade processes in the standard deviation 
of the quantity upon analysis, we reinterpret those mod- 
els considering a cascade processes for the average value 
of the observable, which in this case is the temperature 
anomaly, so that the average value over a certain scale £ 
comes from the multiplicative process, 



h (*)= n^(*^* +i ) 



(18) 



with (II representing the average over a reference period 
L, e.g. L = 11 years, and / representing the fraction 
of measure (average) passing from a earlier generation 



(time scale) to the subsequent. For this case the greater 
i, the smaller £, i.e., i = corresponds to £ = L and £ — n 
yields the scale of the temperature anomaly we are willing 
to describe. Accordingly, the temperature anomaly, £, 
would be equal to, 



& (t) = e (t) + fit (t) , 



(19) 



where e represents a Gaussian independent and identi- 
cally distributed noise with zero average and appropriate 
standard deviation. A skew distribution of /Lt£will lead to 
a skew distribution of £ as well. Concomitantly, in order 
to represent the slight kurtosis we can consider a fur- 
ther contribution, Q, coming from a multiplicative noise 
process, 



C(t) = r, (t) a (t) h (t) 



(20) 



with a (t) being a function of past values of £ and C in 
a heteroscedastic way (see for details and UK) and 
r\ another Gaussian noise not correlated with e. The 
C contribution is modulated by a step function h (t) — 
6 [t - (T - AT)] 6 [T + AT — t] where the center of the 
intervals come from a shot noise following a frequency 
related to Fig. 10. Mathematically, the probability of 
having a contribution from £ arising from a perturbation 
centered at time T is given by, 



p(T)ocy;.*(T-To 



(21) 



and the total temperature anomaly will thus correspond 
to the sum of & (i) and £ (t). 

Alternative models can be presented, namely the mim- 
icking of the fluctuations of the temperature anomaly in- 
stead of that quantity. However, it is simple to find one 
quantity after the other bearing in mind their relation 
and statistical properties as in happens in other problems 
such as fluid turbulence and price dynamics in financial 
markets. To conclude, we would like to mention that the 
mechanisms for the definition of surrogates that we de- 
scribed, namely the phase preservation, can be useful in 
the enquiry into the Gaussian nature of the primary noise 
in processes as those introduced in previously published 
work. 
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